Accurate Effective-One-Body waveforms of 
inspiralling and coalescing black-hole binaries 
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The Effective-One-Body (EOB) formalism contains several flexibility parameters, notably as, 
iipoic and orr. We show here how to jointly constrain the values of these parameters by simulta- 
neously best-fitting the EOB waveform to two, independent, numerical relativity (NR) simulations 
of inspiralling and/or coalescing binary black hole systems: published Caltech-Cornell inspiral data 
(considered for gravitational wave frequencies Muj < 0.1) on one side, and newly computed coa- 
lescence data on the other side. The resulting, approximately unique, "best-fit" EOB waveform is 
then shown to exhibit excellent agreement with NR coalescence data for several mass ratios. The 
dephasing between this best-fit EOB waveform and published Caltech-Cornell inspiral data is found 
to vary between —0.0014 and +0.0008 radians over a time span of ~ 2464M up to gravitational wave 
frequency Muj = 0.1, and between +0.0013 and —0.0185 over a time span of 96M after Mu — 0.1 
up to Muj — 0.1565. The dephasings between EOB and the new coalescence data are found to be 
smaller than: (i) ±0.025 radians over a time span of 730M (11 cycles) up to merger, in the equal 
mass case, and (ii) ±0.05 radians over a time span of about 950M (17 cycles) up to merger in the 
2:1 mass-ratio case. These new results corroborate the aptitude of the EOB formalism to provide 
accurate representations of general relativistic waveforms, which are needed by currently operating 
gravitational wave detectors. 

PACS numbers: 04.25.Nx, 04.30.-w, 04.30.Db 



I. INTRODUCTION 



The Effective-One-Body (EOB) formalism [1-4] is an 
analytical approach which aims at accurately describ- 
ing both the motion of, and the radiation from, coa- 
lescing binary black holes. The EOB method uses, as 
basic input, high-order post-Newtonian (PN) expanded 
results (see [5] for a review of the PN-theory of gravi- 
tationally radiating systems). However, one of the key 
ideas in the EOB method is to avoid using PN results 
in their original 'Taylor-expanded' form (symbolically 
/ Taylor (w/c) = c + Cl v/c+C2{v/c) 2 + - ■ ■ + c n (v/c) n ), but, 
instead, to 're-package' them in some resummed form, 
i.e., symbolically, to replace / Taylor (w/c) by some non- 
polynomial function f E (v/c), defined so as to incor- 
porate some of the expected non-perturbative features 
of the (unknown) result. This re-packaging is crucial 
for being able to bypass the strong limitations of PN 
results. Indeed, by itself PN theory is unable to go 
beyond the (adiabatic) early inspiralling stage of black 
hole coalescence, 1 while the EOB method is able to de- 
scribe, in a continued manner, the full coalescence pro- 
cess: adiabatic early inspiral, nonadiabatic late inspiral, 
plunge, merger and ring-down. The EOB method corn- 



See Appendix B for a new confirmation of this fact 



prises three, rather separate, parts: 

1. a description of the conservative (Hamiltonian) piece 
of the dynamics of two black holes; 

2. an expression for the radiation-reaction force T v that 
supplements the Hamiltonian dynamics; 

3. a description of the gravitational wave (GW) signal 
emitted by a coalescing binary system. 

For each one of these parts, the EOB method uses 
special resummation techniques, inspired by specific re- 
sults going beyond perturbation theory. For instance, 
the resummation of the EOB Hamiltonian (part 1.) was 
inspired by a specific resummation of ladder diagrams 
used to describe positronium energy states in Quantum 
Electrodynamics [6]. The resummation of the radiation 
reaction force was inspired by the Pade resummation 
of the flux function introduced in Ref. [7]. As for part 
3., i.e. the EOB description of the gravitational radi- 
ation emitted by a coalescing black hole binary, it was 
mainly inspired by the classic work of Davis, Ruffini and 
Tiomno [8] , which discovered the transition between the 
plunge signal and a ringing tail when a particle falls into 
a Schwarzschild black hole. 

Before the availability of reliable numerical simula- 
tions, the EOB method made several quantitative and 
qualitative predictions concerning the dynamics of the 
coalescence, and the corresponding GW radiation, no- 
tably: (i) a blurred transition from inspiral to a 'plunge' 
that is just a smooth continuation of the inspiral, (ii) a 
sharp transition, around the merger of the black holes, 
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between a continued inspiral and the ring-down signal, 
and (iii) estimates of the radiated energy, and of the 
spin of the final black hole (the latter estimates were 
made both for nonspinning binaries [2] and for spinning 
ones [9] ) . Those predictions have been broadly confirmed 
by the results of recent numerical simulations performed 
by several independent groups (for a review of numeri- 
cal relativity results see [10]). The recent breakthroughs 
in numerical relativity (NR) [11-15] open the possibil- 
ity of acquiring some knowledge about binary black hole 
coalescence that goes beyond what either PN theory, or 
its resummed avatars (such as the EOB), can tell us. 
Actually, it was emphasized early on [4, 16, 17] that the 
EOB method should be considered as a flexible structure, 
containing several parameters representing (yet) uncalcu- 
latcd results, that would need NR results (or real observa- 
tional data!) to be determined. For instance, Refs. [4, 16] 
introduced a parameter (here denoted as 05) represent- 
ing uncalculated 4 PN, and higher, contributions to the 
crucial EOB "radial potential" A(R). Rcf. [17] intro- 
duced several more EOB flexibility parameters, notably 
Wpoic (entering the Pade resummation of the radiation 
reaction force) and a parameter (here replaced by 5rr) 
describing uncalculated non quasi-circular (NQC) con- 
tributions to the radiation reaction. Recently, Ref. [18] 
augmented the list of EOB flexibility parameters by in- 
troducing two parameters (here denoted as a and b) rep- 
resenting NQC contributions to the waveform, as well as 
two parameters, t m and 5 (together with the choice of 
an integer p), describing the "comb" used in matching 
the inspiralling and plunging waveform to the ring-down 
one. Each one of these EOB flexibility parameters (a 5 , 
^poic, Srr, a, b, t m , 5, p) parametrizes a deformation 2 of 
the originally defined EOB. Each direction of deforma- 
tion, e.g., d/da 5 , hopefully adds some "missing physics" 
that either has not yet been calculated because of tech- 
nical difficulties 3 , or represent only an effective descrip- 
tion of a complicated, nonperturbative process which is 
not directly formalizablc in a calculable way. In both 
cases, the EOB programme aims at using NR results to 
determine the "best fit" values of the flexibility param- 
eters; i.e., the values that, hopefully, allow an analytical 
EOB waveform to accurately represent the exact general 



2 We use here the word deformation in the mathematical sense. 
Ideally we would like the list of EOB flexibility parameters to 
describe a kind of versal deformation of the original EOB, i.e. 
a multi-parameter family which is general enough to encompass 
all the physics contained in real GW coalescence signals, starting 
from the originally defined EOB waveform, which was based on 
a rather coarse representation of the coalescence waveform. 

3 For instance, the exact, general relativistic value of as (or, rather, 
of the ^-dependent coefficient a§{y) = uar, + u 2 a' 5 + ■■■ of 
(GM/c 2 RY' in A(R)) has not yet been calculated simply be- 
cause it would represent a huge technical challenge, involving a 
4 PN (and 4-loop) generalization of the rather involved 3 PN 
(and 3-loop) work that led to the unique determination of the 
lower-order coefficient 04(1/) [5, 19]. 



relativistic inspiralling and coalescing waveform. Note 
that, in this paper, we will not use the terminology of 
faithful (versus effectual) waveforms [7]. Indeed, this ter- 
minology refers to particular measures of the closeness 
of two waveforms (called "faithfulness", J 7 , and "effec- 
tualness", £ in [17]) which are based on specific ways 
of maximizing normalized overlaps. These measures are 
not the best suited for our present purpose because they 
are detector dependent (through the use of the detector's 
spectral noise curve Sh{f) in the Wiener scalar product 
< X,Y >, see e.g., Eqs. (6.1) and (6.2) of [17]). By con- 
trast, we are interested here in hopefully showing that 
EOB waveforms can be "close" to general relativistic ones 
in a much stronger mathematical sense, say in the time- 
domain norm: sup te[tlit2] \ h EOB (t) - h Exact (t)\ < e. 
Actually, the most important "closeness" , for data analy- 
sis purposes, is the closeness of the phases. Therefore we 
shall primarily consider the time-domain phase norm: 
ee inf T , Q su Pte[tljt2] |^ 2 OB (i + x) + « - 0^ 2 R (t) I , 
where we minimize over the two arbitrary parameters r 
(time-shift) and a (phase-shift). When HA^Hoo is smaller 
than e for most physically relevant intervals [£1^2], we 
shall say that the (time-domain) EOB waveform hP OB (t) 
is an e — accurate representation of /i Exact (i). 

The programme of determining the "best fit" flexibility 
parameters by comparing EOB predictions to NR results 
has been initiated in several works [16, 18, 20-22] (see 
also [23-25] for other works involving the comparison of 
EOB waveforms to NR ones). For some parameters, 4 it 
has already been possible to determine them, or, at least, 
to find a rationale that allows one to fix them in a near- 
optimal manner. For instance, it was found in Ref. [18] 
that p = 2, i.e. the use of 2p+ 1 = 5 matching points and 
5 corresponding positive-frequency QNMs was optimal 
from a practical point of view, in the sense that smaller 
values led to visibly worse fits, while higher values led to 
only a rather marginal improvement. We shall therefore 
fix p to the value p = 2. Concerning the "central match- 
ing time" t m , previous work [18, 20, 22] has found that 
it was near optimal to choose (as advocated in [2]) t m to 
be the so-called "EOB light-ring crossing" time, i.e. the 
EOB dynamical time when the EOB orbital frequency fl 
reaches its maximum. Concerning the matching-comb 
width parameter 5 — At/(2p) (where At is the total 
width of the matching interval), Refs. [18, 22] found that 



4 Note that several of the EOB flexibility parameters listed above 
refer to the recently introduced resummed 3 +2 -PN accurate 
EOB waveform [18, 21] and to the "comb" matching proce- 
dure of Rcf. [18]. The EOB dynamics and waveform used in 
the works of Buonanno and collaborators differ in several ways 
from the dynamics and waveform used by us, notably: (i) a ra- 
diation reaction force of lesser PN accuracy, and without NQC 
corrections, (ii) a waveform of "Newtonian" accuracy without 
NQC corrections, and (iii) a simpler matching procedure between 
the plunge and the ring-down involving only three (positive- 
frequency) quasi-normal modes (QNM) and an instantaneous 
matching (as used in some of the original EOB papers [2, 26]). 
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S = 1.7 Mfi na i (corresponding to At = 45 = 6.8M final ) 
yielded a good result. Here Af nna i denotes the mass of 
the final black hole. Here also, we fix t m = iEOjT"" 8 ' 
and 5 = 1.7Af nna i. Moreover, we shall discuss below a 
rationale allowing one to fix the parameters a and b. 

Summarizing: the only EOB flexibility parameters 
which have not yet been uniquely determined are a^, 
fpoio and cirr. Some recent works [20-22] have explored 
the influence of these parameters on the EOB waveform 
and have made initial steps towards determining 'best 
fit' values for these parameters; i.e., values leading to an 
optimal agreement with NR data. In particular, Ref. [20] 
found that the faithfulness T (in the sense of Sec. VIA 
of Ref. [17]) of restricted EOB waveforms against NASA- 
Goddard NR coalescence waveforms was largest when 5 
as belongs to some rather wide interval, say [20, 100], cen- 
tered around as ~ 60. On the other hand, Ref. [21] found 
that the accuracy (in the sense of the L x norm of the 
phase difference) of the resummed 3 +2 -PN EOB wave- 
form 6 with respect to the high-accuracy Caltech-Cornell 
(CC) NR long-inspiral waveform was at its best when 
as belonged to an interval ~ [10,80] centered around 
a 5 ~ 40. The influence of the flexibility parameter w po io 
was studied in Refs. [21, 22]. It was found that, for 
any given values of 05 and 5rr, and for any given NR 
waveform, there existed a well determined value of v po \ e 
that minimized the phase difference between EOB and 
NR (see below for a more precise formulation). How- 
ever, those previous EOB-NR comparisons limited them- 
selves to considering one NR data set at a time (the 
published Caltech-Cornell inspiral data for Ref. [21], and 
some Albert Einstein Institute (AEI) coalescence data 
for Ref. [22]). 

The aim of the present paper is to go beyond this 
piece-meal consideration of NR data and to best fit (in 
phase) the flexed EOB waveform, /i EOB (as, w po ic, orr; t), 
simultaneously to several independent NR waveform data 
(namely inspiral and coalescence data produced by the 
Jena group and reported here, and published inspiral 
Caltech-Cornell data). Our main result will be that the 
best fit values of the three remaining EOB flexibility pa- 
rameters (a5,i> po io, orr) are approximately determined, 
in the sense that they must all take values in relatively 
small, correlated, intervals. It is then found that the 
resulting, approximately unique, best fitted EOB wave- 
form exhibits a remarkable agreement (modulo differ- 
ences compatible with estimated numerical errors), both 
in phase and in modulus, not only with the data that we 
use in the fit (i.e., equal-mass Jena data and equal-mass 
Caltech-Cornell data considered for Muj < 0.1), but also 
with other NR data (namely, unequal-mass Jena data and 



5 Note that Ref. [20] uses the notation A for 05. 

6 We refer to the PN accuracy of this waveform as 3 +2 PN because 
it includes not only the known comparable mass 3 PN waveform 
corrections, but also the test-mass limit of the 4 PN and 5 PN 
waveform amplitude corrections [21]. 



Caltech-Cornell data after Muj = 0.1). 

Our work focusses on the comparison between the EOB 
predictions and NR data because the EOB method is 
the only existing analytical approach which: (i) incorpo- 
rates, in an exact manner, all the theoretical knowledge 
acquired through many years of post-Newtonian stud- 
ies, (ii) provides waveforms covering the full coalescence 
process from early inspiral to ring-down, and (iii) can de- 
scribe spinning binaries (see, in this respect Refs. [4, 27]). 
However, as some studies have emphasized the nice prop- 
erties of one specific PN approximant, called TaylorT4 
in [28] (for consistency with the Tl, T2 and T3 Tay- 
lor approximants considered in [29]), we shall discuss it 
briefly in Appendix B, though it docs not satisfy our 
requirements (ii) above, namely that of providing wave- 
forms covering the full coalescence process. 

This paper is organized as follows. In Sec. II we briefly 
describe the numerical simulations, whose results we use 
in the following. Section III summarizes the definition of 
the 3 +2 -PN accurate EOB waveform that we use. Sec- 
tion IV is the central section of this work: it shows how 
the simultaneous comparison of EOB to two different 
NR data sets allows one to determine a small range of 
'best fit' (correlated) EOB flexibility parameters 05, v po i c 
and (Zrr. Section V selects central values for the best 
fit parameters and discusses in detail the agreement be- 
tween the EOB waveform and the Jena NR waveform, 
for various mass ratios. The paper ends with a conclud- 
ing Section, followed by two Appendices. Appendix A 
is devoted to the issue of determining the metric wave- 
form h(t) from the curvature waveform ipi(t), while Ap- 
pendix B discusses the TaylorT4 approximant. Except 
when otherwise specified, we use in this paper units such 
that G = c = 1. 



II. NUMERICAL RELATIVITY SIMULATIONS 

Numerical simulations were performed with the BAM 
code [30, 31], which evolves black- hole binaries using 
the "moving-puncture" approach [12, 13]. The relevant 
physical and numerical parameters for our simulations 
are presented in Table I; note that the results from the 
equal-mass simulations were presented in [32] , which also 
contains extensive error analysis and comparison with 
standard post-Newtonian inspiral approximants. These 
results are also in good agreement with those of [28] 
over the shared frequency range. We shall present be- 
low an explicit comparison of the phase of the waveform 
of Ref. [28] with the one of our equal-mass simulation. 

A. Initial data 

Following the moving puncture approach we set up ini- 
tial data containing two black holes via a Brill-Lindquist- 
like wormhole construction [33], where the additional 
asymptotically flat end of each wormhole is compacti- 
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fied to a point, or "puncture". The entire 2-wormholc 
topology can thus conveniently be represented on Br. It 
has long been understood how to set up such puncture 
initial data, and in particular how to avoid working with 
divergent quantities [34-37] . More recently it has turned 
out that the gauge conditions used in the moving punc- 
ture approach actually allow a simpler representation of 
the black hole interior during the evolution: the black- 
hole throat is pushed an infinite proper distance away 
from the horizon, and the initial puncture geometry is re- 
placed by a new compactified asymptotics with a milder 
singularity [38-40]. 

One key element of the simplicity of the moving punc- 
ture approach is that black holes can be modeled on a 
Cartesian numerical grid without the need to deal with 
black hole excision techniques. Another is that the as- 
sumption of an initially conformally flat spatial geometry 
yields a very simple way to generate any number of mov- 
ing, spinning black holes [36, 41]. Note however that the 
puncture initial data are not restricted to conformal flat- 
ness a priori [34, 35, 37], and generalizations that better 
model spinning black holes have been suggested [42, 43]. 

Assuming conformal flatness for the initial data, and 
assuming the extrinsic curvature of the initial slice to 
be within the class of nonspinning Bowen-York solu- 
tions, the freedom in specifying initial data comprises 
the masses, locations and momenta of each black hole. 

The mass of each black hole, Mi (i — 1,2), is speci- 
fied in terms of the Arnowitt-Deser-Misner (ADM) mass 
at each puncture, which is, to a very good approxima- 
tion [44-46] equal to the irreducible mass [47, 48] of the 
apparent horizon 

where A$ is the area of the apparent horizon. We iden- 
tify this mass with the mass (denoted below mj) used in 
post-Newtonian theory. This assumption will only hold 
exactly in the limit where the black holes are infinitely 
far apart and stationary, but we consider any error in 
this assumption as part of the error due to starting the 
simulation at a finite separation. 

The constraint equations for black-hole binary punc- 
ture initial data are solved using a pseudo-spectral code 
[49], and resampled for our finite difference grid using 
high-order polynomial interpolation [31]. 

We want to specify initial data for non-spinning black 
holes in the center-of-mass frame, such that the trajec- 
tories correspond to quasicircular inspiral, i.e. the mo- 
tion is circular at infinite separation, and the eccentric- 
ity vanishes. Following [50], we obtain the initial mo- 
menta of the black holes from a post-Newtonian inspiral 
calculation, using a 3PN-accurate conservative Hamilto- 
nian [19], and 3.5PN accurate beyond leading order orbit- 
averaged radiation flux [51, 52]. We have measured the 
eccentricity from oscillations in the separation and fre- 
quency as described in [50], and have obtained the val- 
ues 0.002, 0.003, 0.005, for mass ratios q = 1,2 and 4 



respectively. 

B. Numerical evolution 

We use the BSSN formulation of the Einstein equations 
[53, 54] for time evolution, which are formulated in terms 
of a conformal 3-metric 7^ , related to the physical metric 
as 

lij=Xlir ( 2 ) 

Representing the conformal factor by the quantity \ nas 
the advantage that, when dealing with puncture data, 
the conformal factor \ conveniently vanishes at each 
puncture [12]. Details of our implementation of the 
BSSN/moving-puncture system are described in [30]. We 
also need to choose a lapse and shift during the evolution 
to determine our coordinate gauge. As is common in the 
moving puncture approach, we use the "1+log" slicing 
condition [55] 

d a = -2aK, (3) 
and the f-driver condition [56, 57], 

d P l = \b\ (4) 
d Q B l = d r- V B\ (5) 

where do — d t — fi l di. The parameter 77 in the shift- 
condition effectively regulates the coordinate size of the 
apparent horizons, and is set to r\ = 21 M in our simula- 
tions. 

The Einstein evolution equations are solved numer- 
ically with standard finite-difference techniques as de- 
scribed in [30, 31]. Spatial derivatives are approximated 
with sixth-order accurate stencils. First order derivatives 
corresponding to Lie derivatives with respect to the shift 
vector are approximated with off-centered operators as 
described in [31], all other derivatives are approximated 
with centered finite difference operators. Kreiss-Oliger 
artificial dissipation operators which converge to zero at 
fifth order are applied as described in [30, 31]. Time 
evolution is performed with a fourth-order Runge-Kutta 
integration. Our box-based mesh refinement is described 
in [30] . Time interpolation errors in the mesh-refinement 
algorithm converge only at second order, but do not seem 
to contribute significantly to the error budget, as does the 
Runge-Kutta time integration. 

The grid configurations we have used for our equal 
mass runs are described in [32]. For the unequal mass 
runs, we have used the 56,64,72-gridpoints configurations 
of [32], adding two further refinement levels to push the 
outer boundary further out by roughly a factor of four. 
In order to be able to re-use our equal mass grid config- 
urations, we always choose the mass of the smaller black 
hole, which determines our resolution requirements, at 
Mi = 0.5. 
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C. Wave extraction 

The gravitational wave signal is extracted at differ- 
ent surfaces of constant radial coordinate by means of 
the Newman-Penrose Weyl tensor component ^4 [58, 59] 
which is a measure of the outgoing transverse gravita- 
tional radiation in an asymptotically flat spacctime. At 
finite distance to the source the result depends on the co- 
ordinate gauge and the choice of a null tetrad. Our choice 
of tetrad and details of the wave extraction algorithm are 
described in detail in [30] . We choose our extraction sur- 
faces at 40, 50, 60, 80, and 90 M. In [32] we extrapolated 
the waveform amplitude (though not its phase) to the 
value that would be observed at infinity; in this work we 
deal with the raw numerical data at the farmost extrac- 
tion radius, but use some extraction-radius-extrapolatcd 
results to provide uncertainty estimates. See next sub- 
section. 

The analysis carried out in this paper will use, as ap- 
proximate asymptotic waveform, the curvature perturba- 
tion extracted at radius 90M, without any extrapolation 
(neither with respect to extraction radius, nor with re- 
spect to resolution) . The comparisons between numerical 
data and analytical predictions done below will make use 
of metric (by contrast to curvature) waveforms. We dis- 
cuss in Appendix A the integration procedure we used 
to compute the numerical metric waveform from the raw 
curvature waveform output of the simulations. In this 
paper, we focus on the i = m = 2 "quadrupolar" wave- 
form. 



D. Accuracy 

1. The equal-mass case 

A detailed error analysis was performed for the equal- 
mass waveforms in [32] . In this section, we will first sum- 
marize the results of that error analysis, and then com- 
plete it by more carefully quantifying the uncertainty in 
the phase. As we shall see, our refined estimate of the 
uncertainty in the phase will end up being significantly 
lower than the upper bound, 0.25 radians, quoted in [32]. 

The amplitude and phase of the equal-mass nonspin- 
ning waveforms show sixth-order convergence with re- 
spect to numerical grid resolution prior to merger, with 
a small drop in convergence order around merger time. 
Higher-accuracy results were constructed by Richardson 
extrapolation with respect to numerical resolution, and 
this procedure also allowed an estimation of the contri- 
bution to the uncertainty in the amplitude and phase 
from discretization error. The discretization error in the 
amplitude was found to be below 0.5%, while the dis- 
cretization error in the phase was estimated to be below 
0.01 radians. These are conservative error estimates ob- 
tained by observing the numerical errors over the course 
of the entire simulation. See in particular Fig. 4 in [32]. 
If we look at that figure we may conclude that the error 



estimate of the phase is extremely conservative, but one 
should also be aware that the quoted numerical phase 
error takes into account only instantaneous differences in 
the value of the waveform phase, but not secular drifts. 
When the analysis for [32] was performed, the authors 
hoped that the conservative value quoted would account 
for any phase drifts. We shall see below that, however, 
there might remain sources of secular drifts that are not 
yet well understood. 

In addition to the discretization error, there is also an 
error due to measuring the waveform at a finite distance 
from the source. For both the waveform amplitude and 
phase, it was found in [32] that finite extraction radii er- 
rors were much larger than discretization errors. Prior 
to merger, the error in the amplitude was found to fall 
off as l/-Rg X , where R ex was the radiation extraction ra- 
dius, and this observation allowed a clean extrapolation 
to R ex — > 00, and, once again, an estimate of the uncer- 
tainty in the amplitude. The uncertainty in the extrapo- 
lated amplitude was at most 2% before merger. Around 
merger time, the amplitude error fall-off is dominated by 
a 1 /R ex term, and the uncertainty in the extrapolated 
amplitude grows to around 5%. However, in this paper 
we use the raw data calculated at the extraction radius 
R ex = 90M, and as such the uncertainties are larger, as 
much as 5% over the entire simulation. The largest un- 
certainties in the finite-extraction-radius amplitude are 
at early times, when the amplitude is small, and around 
merger, when the dynamics are strongest. 

In [32] the total phase uncertainty accumulated on a 
time interval of duration 1400M extending up to gravi- 
tational wave frequency Muj = 0.1 was quoted as being 
0.25 radians. This large value was an upper bound which 
was quoted in view of the difficulty in finding a robust 
method to extrapolate the phase to infinite extraction ra- 
dius. These difficulties were related to the specific phase 
alignment method which was used in [32]. There, one 
was first choosing some frequency at which to line up the 
phases and frequencies of waves from different extraction 
radii, and then attempting to perform an extrapolation. 
Although it is entirely valid to time- and phase-shift any 
number of waveforms to perform a comparison between 
them, it turned out that this is not an efficient way to 
perform a consistent extraction-radius extrapolation. 

By contrast, for the purpose of the present paper we 
have performed a new study of the extraction-radius ex- 
trapolation which follows the strategy proposed in [28]. 
More precisely we used two similar, but different, phase 
alignment methods. The first one consists of simply in- 
troducing the "Newtonian retarded time", at the coor- 
dinate extraction radii R ex , un = t — R ex and study 
the waves as function of %. Then, when attempt- 
ing extrapolation with respect to R ex , we find a clear 
Cq(un) + C2(wn) / R\ x fall-off in the error, and are able to 
make a clean extrapolation to infinity. We have also re- 
peated the analysis with an extra +C3(un) / 'R\ x term in 
the fit. 

The second method consists of using, inspired by the 



6 



TABLE I: Details of the simulations discussed in the test. From left to right, the columns report: mass ratio q — 777,2/777,1; 
symmetric mass ratio v — m\m2/(m\ + m'2) 2 ; initial coordinate separation D of the punctures; inital ADM mass; initial 
tangential (pt) and radial (p r ) momentum of the black holes; mass and dimensionless spin parameter jf = Jt/Mf of the 
final black hole; mass and dimensionless spin parameter J™ 5 and M ( rmg /M of the final black hole obtained only from the 
ringdown; dominant (quasi-normal-mode) frequency of the ringdown. Quantities are scaled by the total initial black hole mass 
M = mi + 7712 as indicated. 
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result in [28], the a priori more accurate definition of 
retarded time, uq — t — r*, where the (approximate) 
Regge- Wheeler tortoise coordinate r* is (following [28]) 
defined as = R ex + Madm + 2Madm log[(R ex + 
-Madm)/(2Madm) — 1] ■ This improved choice of retarded 
time allows us again to perform a clean extrapolation 
to infinity. As when using itjvj, we use two different 
fits: c (u B ) + c 2 (u B )/Rl x and c (u B ) + c 2 (u B )/R 2 ex + 
c 3{ u b)/ Rex- We then estimate the uncertainty in the 
phase of the farmost unextrapolated data, extracted at 
Rex — 90M, by comparing the following five phases: (i) 
the raw phase 4>qqm measured at R ex = 90M, (ii) the 
phase </>^2 extrapolated using «n and assuming a \jR\ x 
fall off, (iii) the phase 4>^ 3 extrapolated using «n and 
assuming a 1/R% X + l/R^x ^ au on? ! tne phase (j) B2 
extrapolated using ub and assuming a 1/ R% x fall off and 
(v) the phase </>^? 3 extrapolated using hb and assuming 
a l/Rg X + 1/Rlx fall off. The differences between the 
phases are computed after they have been aligned by us- 
ing the two-times pinching technique of Ref. [22] (which 
is reviewed in Sec. IV below). For consistency with our 
EOB-NR matching discussed in Sec. V below we use 
as "pinching" gravitational wave frequencies wi « 0.1 
and u>2 = 0.4717. Note that these frequencies bracket 
the merger time. The four phase differences 0goM — 4>T 
where i £ {N2, N3, B2, B3} are exhibited as functions 
of the numerical relativity coordinate time at 90M, in 
Fig. 1 . The triangles in the figure indicate the two times 
corresponding to the two pinching frequencies (ui,^), 
while the vertical dashed lines indicate the time interval 
[1200, 1900] w [t;l, £r] which will turn out to be crucial for 
our analysis in Sec. V below. Several conclusions can be 
drawn from Fig. 1: First, the choice of retarded time, «n 
or ub, does not matter much for the extrapolation pro- 
cedure. Second, though the phase differences over the 
entire span of the simulation can reach values ~ +0.13 
radians around merger time (peak at t ~ 1930M) and/or 
~ —0.2 radians (during ringdown), they stay quite small 
during the time interval [ti,, £r] that we shall focus on in 
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FIG. 1: Differences between the phase extracted at R £X = 
90M, and the phase extrapolated to infinity based on two 
choices of the retarded time and on two choices of the extrap- 
olating polynomial, as described in the text. The choice of 
retarded time makes little difference to the result. 



our analysis below. Most importantly for the following 
the maximum phase differences over the interval [t;l,^r] 
stay within the rather small interval (—0.042, +0.032) ra- 
dians. 



7 Note also that with the above choice of pinching times (wi, W2) 
the phase differences stay quite small, namely —0.06 radians, 
during the entire inspiral. However, this result sensitively de- 
pends on the way the phases have been matched. For instance, 
when using pinching frequencies w\ = u>2 = 0.1 one observes 
maximum phase differences of ~ +0.07 radians at merger and 
~ —0.45 radians during ringdown, while they stay between (- 
0.03,0) radians during the inspiral. On the other hand, when 
using pinching frequencies around merger, i.e., u>\ = 0.36 and 
OJ2 = 0.38, one gets quite small phase differences during merger 
and ringdown, but one observes large dephasings at early times, 
that build up to —0.6 radians. 
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FIG. 2: Comparison between Caltech-Cornell and Jena actual 
numerical data: the phase difference A022 CJena = <t>° C - </>9om 
is shown versus Caltech-Cornell GW frequency tvcc- 



After this "internal" way of trying to estimate the nu- 
merical errors in the phase of our equal-mass simulation, 
let us consider an "external" way which consists in di- 
rectly comparing the unextrapolated, 90M raw phase 
</>90M (i) to the phase computed by Boyle et al. [28] and 
kindly communicated to us. In Fig. 2 we are directly com- 
paring two phases: our unextrapolated (j>goM(t) and the 
resolution- and radius-extrapolated Caltech-Cornell cur- 
vature phase 4>cc(t)- The phase difference A02 2 CJcna = 
tftcc — 09OM was plotted versus the Caltech-Cornell (cur- 
vature) frequency ujcc- This phase difference was com- 
puted in the following way. First, we used the two- 
pinching frequencies u>i = 0.059517 and 0J2 — 0.14976 
(indicated by two dashed vertical lines in the figure) 
to determine the time and phase shifts (r, a), see be- 
low, then the result A^ 2 CJcna (t cc ) = ^cc(icc) - 
(<^90M(icc + t) + a) is plotted versus wcc instead of 
tec- In addition, since the Caltech-Cornell simulation 
extends only up to usqq* ~ 0J2 ~ 0.15, we have estimated 
three different possible extrapolations of the phase dif- 
ference A<^ 2 CJcna beyond u 2 and up to lu r = 0.1898. 
8 These three different extrapolations were obtained by 
fitting A</>(wcc) over the interval [0.1,0.15] by three dif- 
ferent polynomial functions of wcc : quadratic, cubic and 
quartic. As we see on Fig. 2, the quadratic fit is the one 
which gives the worst possible phase difference over the 
interval [o>l,wr,]. We use this worst case as estimate of 
the maximum phase difference between Caltech-Cornell 



8 Note that 0Ji = ojl corresponds to the lower limit of the 
crucial EOB-NR comparison interval used in Sec. IV, while 
corresponds to its upper limit tj^. 



and Jena phasings over [fi,, £r] ■ More precisely, while the 
minimum value of A(/>2 2 CJona over the interval is —0.0068 
radians, its maximum is +0.04484 radians (at wr, i.e., 
at the extreme right of Fig. 2). This corresponds to a 
two-sided CC-Jena phase difference (in the sense of foot- 
note 12 of Ref. [22] ) ±l/2(0.04484-(-0.0068)) = ±0.026 
radians over the interval [wl, wr]. As this difference a pri- 
ori comprises many possible "noise sources" coming from 
comparing two very different simulations, with different 
wave extraction procedures, we consider that this is our 
best present estimate of the unknown "real" error on the 
difference (j>goM — 4>cc- In addition, as a recently pub- 
lished upper limit on the total accumulated phase error in 
the Caltech-Cornell data of Ref. [28] is 0.01 radians [60] 
(which is rather small), we shall consider in the following 
that ±0.026 radians provides our best current estimate 
of the real error on the equal mass 4>qom over the in- 
terval [wl,wr]. Note, in passing, that the internal error 
analysis procedure discussed above was giving a roughly 
comparable error estimate, namely a two-sided phase dif- 
ference ±1/2(0.032- (-0.042)) ~ ±0.037 radians. How- 
ever, we cannot rely on this internal analysis because it 
fails to explain the origin of a striking feature of Fig. 2, 
which is that, before the plateau of very small phase dif- 
ferences reached between frequencies 0.06 and 0.15, there 
is a steeper phase gradient which reaches —0.16 radians 
at oj = 0.04445, roughly corresponding to the beginning 
of the Jena simulation. Part of this error may be due 
to residual eccentricity — a quick comparison with post- 
Newtonian results using the techniques described in [50] 
suggests that the phase error from a residual eccentricity 
of e ~ 0.002 could be as much as 0.05 radians. We feel, 
however, that most of the error comes from some secular 
drift at early times which is not yet well understood. 



2. The unequal-mass cases 

For the unequal-mass case 2:1, we find similar results, 
namely, that the finite extraction radii dominate the er- 
ror, and the amplitude error is below 5% prior to merger. 
As for the accumulated phase error in the 700M time 
span up to Mlu = 0.1, we did not carry out the ra- 
dius extrapolation analysis discussed above in the un- 
equal mass case. As a rough upper limit we quote an 
accumulated phase error of 0.15 radians. In contrast to 
the equal-mass case, the fall-off in the amplitude error 
with respect to radiation extraction radius is not so clean 
around merger time, preventing us from performing an 
accurate extrapolation to infinity. As such, we would 
conservatively give an uncertainty estimate of 10% of the 
amplitude at merger and later. 

In the unequal-mass case 4:1, the case is different 
again: here the discretization error dominates the phase 
error, suggesting that higher-resolution simulations are 
needed. Our estimate for the accumulated phase uncer- 
tainty up to Mui = 0.1 is 0.25 radians, based entirely 
on discretization error. For the amplitude we estimate 
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that the uncertainty is similar to that in the 1:2 case, 
i.e., around 5% prior to merger, and 10% after that time. 



E. Final parameters of the black hole 

The final mass of the black hole is obtained by sub- 
tracting the radiated energy from the initial mass. While 
the initial mass (the ADM mass) is known very accu- 
rately from the solution of the constraints with spectral 
methods [49], the radiated energy is less accurate and 
dominates the errors of the final mass and Kerr spin pa- 
rameter. The radiated energy is not very accurate, due 
to the loss of accuracy in the wave signal at merger time 
for the equal mass case (leading to a conservative error 
estimate of 4%), and the problems of extrapolation in 
radius and gridspacing for the unequal mass cases, which 
lead us to a conservative error estimate of 10% in those 
cases. 

The error in the radiated energy also dominates com- 
puting the quantity jf — J/M 2 , where we either compute 
J from a surface integral as in [30] and M as described 
above, or we calculate jf itself from the ringdown. The er- 
ror in computing the angular momentum J from a surface 
integral falls off very quickly with separation. The dom- 
inant error in this quantity comes from high-frequency 
numerical noise in the integrals, which is however much 
smaller than the error in jf resulting from errors in the 
final mass. 

To determine the mass and spin parameter of the fi- 
nal black hole from the ringdown, we have performed 
two types of fits to the dominant mode. First, the qual- 
ity factor has been obtained from a fit to the domi- 
nant quasi-normal mode 9 of the complex ringdown wave- 
form. This fit was performed by a non-linear least- 
squares Gauss-Newton method, using exp(— at + p) as 
a parameter-dependent template (with two complex pa- 
rameters (a, p)), and an appropriate time interval during 
the ringdown (chosen by minimizing the post-fit resid- 
ual). Then, from the best-fit value of a (i.e., the QNM 
dominant complex frequency <r^ 2 o )> we computed the 
values of (M f rmg , j" ng ) of the final black hole by using 
the interpolating fits given in Appendix E of [61]. The 
triplets (M f ring /M, j f ring ,Mcr+ 20 ) are listed in Table I. 

This method does not require knowledge of the final 
mass, but is actually not well conditioned due to the 
shape of the curve j(co). Better accuracy is obtained by 
only using the real part of the frequency, then again, the 
error in j is dominated by the error in the radiated en- 
ergy. The values are consistent with the values obtained 
from the surface integrals for the angular momentum J. 
The numbers Mf and jf quoted in Table I are consistent 



9 In the notation introduced in Sec. Ill below, the dominant mode 
corresponds to the labels (±,£,£',m,n) = (+,2,2,2,0). 



with both methods, and with the analytical fit for jf pub- 
lished for shorter and less accurate waveforms in [62] . By 
contrast (M f ring , j f ring ), are "best-fit" values that will be 
used below to compute the EOB ringdown waveform. 



III. THE EOB WAVEFORM 

We shall not review here the EOB formalism, which 
has been described in several recent publications [20- 
22, 24, 63]. We refer to these papers, and notably to 
Rcfs. [21, 24], for detailed definitions of the dynamics 
and of the waveform. Let us only indicate here a few 
of the crucial elements of the EOB implementation that 
we use here. We recall below the main ingredients of 
the EOB formalism, focusing on the appearance of the 
various EOB flexibility parameters. 

• The EOB Hamiltonian i? roa i describes the conser- 
vative part of the relative two-body dynamics. We 
use for the crucial "radial potential" A(r) entering 
this Hamiltonian the Pi Pade resummation of 

A Taylor (a 5 , u; u) = 1 - 2u + 2vu z + a±vu A + a b vu b , (6) 

where [3, 19] 

94 41 2 

° 4 = T _ 32 7r ' (7) 

where as is the 4 PN flexing parameter introduced 
in [4], and where 10 u — 1/r. 

• The EOB radiation reaction force ^(wpoiej Srr, v), 
that we shall use here, has the form 

•7>(v P oic, «rr, v) = /^ C (a RR )J c '°(wpoio, v), (8) 

where J^(v po ic, v) is defined MS cl Pade resumma- 
tion [7] of its Taylor expansion. See Eq. (17) 
of [18] where /dis is the P4 Pade resummation of 
(1 - v/v polc )F Ta y lol (v;v). In addition, the factor 
/rr ° 1S a non quasi-circular (NQC) correction fac- 
tor of the form 

This factor was introduced in Refs. [18, 22] (see also 
Ref. [17]). We fix the value of crr to err = 0.2 as 
in [22]. 



Except when said otherwise, we henceforth systematically scale 
dimensionful quantities by means of the total rest mass M = 
mi + ni2 of the binary system. For instance, we use the dimen- 
sionlcss EOB radial coordinate r = R^ob/M, with G = 1. Note 
also that v = n/M with fi = m\milM. 



9 



• We use improved "post-post-circular" EOB dy- 
namical initial data (positions and momenta) as 
in [21, 22]. 

• We use the resummed 3 +2 PN accurate "inspiral- 
plus-plungc" Zcrilli-Moncrief normalized metric 
waveform introduced in Ref. [18, 21]. It has the 
form 

%T lunsc (a,b; u,t) = -4^u(r u nffg QC (a,b)F 22 (u) ( 

(10) 

Here $(f) is the EOB orbital phase, = $ is 
the EOB orbital frequency, r w = r[^(r,p v )] 1 / 3 is 
a modified EOB radius, with being defined as 



-2i# 



dA(r) 
dr 



l + 2v 



Mr) 1 + P J 



(11) 



which generalizes the 2PN-accurate Eq. (22) of 
Ref. [26]. The factor F 22 is a resummed, 3 +2 -PN- 
accurate complex amplitude correction valid during 



the (adiabatic) inspiral (see [21]), and / 22 M (a, b) 
is the following extra complex correcting factor, 
aimed at taking care (in an effective way) of various 
non quasi-circular effects during the plunge 



/ 2 N 2 QC M) 



i 



(rO) 2 



\h Pr * 



(12) 



where p rt is the momentum conjugate to the EOB- 
tortoisc radial coordinate r*. Here we shall fix 
e a = 0.12. In these equations, we have only indi- 
cated the explicit appearance of the waveform flex- 
ibility parameters (a, b). In addition, the waveform 
is, evidently, implicitly depending on a$, which en- 
ters the Hamiltonian, as well as on w po i and arm, 
that enter the radiation reaction force. 



ordinary spherical harmonics Yg m (6, <p) (used as ex- 
pansion basis to define ^ i m ) can be expanded in the 
spheroidal harmonics Sei m (aa,9)e lm ''' characteriz- 
ing the angular dependence of the Kerr-background 
QNMs [64]. As explained in Sec. Ill of Ref. [22], we 
use five positive frequency QNMs computed start- 
ing from the values of M" ng /M and j f nng listed in 
Table I. 

• We match the inspiral-plus-plunge waveform to the 
ring-down one, on a (2p + l)-tooth "comb" (t m — 
pS, t m — (p—l)6, ...,t m — 8, t m ,t m + 5, . . . ,t m +p5), 
of total length At = 2pS, which is centered around 
some "matching" time t m . We fix the integer p 
to the value p — 2, corresponding to five match- 
ing points. As mentioned above, we follow previ- 
ous work [2, 22, 23] in fixing the "matching time" 
t m to coincide with the so-called "EOB light-ring" , 
i.e. the instant when the orbital frequency Cl(t) 
reaches its maximum (this defines, within the EOB 
approach, the merger time). As in [22], we fix 
5 = 1.7M f rmg , which corresponds to a total width 
for the matching interval At = 45 = 6.8M" ng . 

• Finally, we define the complete EOB matched 
waveform (from t = — oo to t = +oo) as 

*22° B ( a 5, "polo, Srr, a, 6, t m , 5; v,t) 
ee 6(t m - t)V™ plmsc (t) + 6(t - t m )^ down (t), (14) 

where 0(t) denotes Heaviside's step function. 

This metric EOB waveform then defines a corre- 
sponding curvature waveform, simply (modulo a 
factor r and normalization conventions) by taking 
two time derivatives of (14), namely 



dt 2 



dt 2 



(15) 



• We use a ringdown waveform, 

*™ gdown (*) = £ c $ e " < * < 



(13) 



N 



where the label N actually refers to a set of in- 
dices (£,£', m, n), with (£,m) — (2,2) being the 
Schwarzschild-background multipolarity degrees of 
the considered <Jv TO waveform with n = 0,1,2,... 
being the "overtone number" of the considered 
Kerr-background Quasi-Normal Mode (QNM; n = 
denoting the fundamental mode), and £' the 
degree of its associated spheroidal harmonics 
Sfm {a<J, 9) . In addition o~~^ = + iuj~^ refers to 
the positive complex QNM frequencies (a^ > 
and uj^ > indicate the inverse damping time 
and the oscillation frequency of each mode respec- 
tively) . The sum over £' comes from the fact that an 



where N e = (£, + 2){£ + !)£{£ - 1) (sec Ap- 
pendix A). Note, however, that in view of the im- 
perfect smoothness 11 of the EOB matched metric 
waveform (14) around t — t m , we find it more con- 
venient, when comparing EOB to numerical data 
that include the merger, to work with the metric 
waveform without taking any further time deriva- 
tives. We discuss in Appendix A the procedure 
that we use to compute from the numerical rela- 
tivity curvature waveform a corresponding metric 
waveform by two time integrations. 



11 A partial cure to this problem would consist in replacing the dis- 
continuous step function 8(t — t m ) in Eq. (14) by one of Laurent 
Schwartz's well-known smoothed step functions (or "partitions 
of unity") 0e((t-t m )/(2p8)). 
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A. Fixing the a and b flexibility parameters 

In this brief subsection we discuss a rationale for choos- 
ing two of the EOB flexibility parameters mentioned 
above, namely a and b, that enter the NQC waveform 
correction factor (12). 

Ref. [18] found that it was near optimal to fix the NQC 
parameter a entering the modulus of the waveform 12 so 
as to ensure that the maximum of the modulus of the 
EOB quadrupolar metric waveform sits on top of that 
of the EOB orbital frequency, i.e., at the "EOB light- 
ring" 13 . We shall therefore "analytically" determine the 
value of the waveform NQC parameter a, as a function 
of the symmetric mass ratio v = p/M = m\mij{m\ + 
77i2 ) 2 by imposing the following requirement: that the 
maximum of |^ , 22 >B (0| be on top of the f2(f). 

In principle, the determination of a by this require- 
ment depends on the choice of the other EOB flexibility 
parameters. In other words, the satisfaction of this con- 
dition will determine a as a function of all the parame- 
ters entering the EOB dynamics and inspiral waveform: 
a = 01(015, <2rr, Vpoiei v)- In practice, however, and as a 
first step towards a fully consistent choice of all the EOB 
flexibility parameters, we fixed a in the following way. 
In previous work it was found both analytically (when 
v <C 1, see Ref. [18]) and numerically (when v = 1/4, 
see Ref. [22]) that the value a = 0.5, together with 
e a = 0.12, led to a sufficiently accurate solution of the 
above requirement. For the present work, we partially 
took into account the parameter dependence of a by fix- 
ing (a$, Vp i e , 5rr) to the central best-fit values that we 
will select below and by then numerically finding the op- 
timal value of a as a function of v only. In particular, we 
identified the following pairs (v, a) of near-optimal val- 
ues: (0.25,0.44), (0.2222,0.49), (0.16,0.64), (0.05,0.905) 
and (0.01,0.985). These are the values that we shall use 
in this work. Note also that the ^-dependence can be 
approximately represented by a simple linear fit, namely 
a{v) = 1.019 - 2.345za 

As for the NQC parameter b entering the phase of the 
(quadrupolar) waveform, previous work [22] has found 
that it had a very small effect (when using the new, 3 +2 - 



12 Here a and b denote the parameters called a' and b' in footnote 
9 of [18] 

13 Note that this coincidence in the locations of the maximum 
of | /i22 (*) I and of fi is automatically ensured when one uses 
(as advocated in [2]) a "restricted" EOB waveform ^/22(t) oc 
Q 2 / 3 exp[— 2i<E>(t)]. It is, however, a non trivial fact that NR 
results show (both in the test-mass limit [18] and in the equal- 
mass case [22]) that the maximum of |*22 (t)| does occur very 
near the maximum of the (corresponding, best matched during 
inspiral) EOB orbital frequency Q(t). This can be considered as 
another successful prediction of the EOB formalism. Note that 
this property does not apply to the maximum of the modulus of 
other GW quantities, such as the instantaneous energy flux or 
the modulus of quadrupole curvature waveform r-0| 2 (t), which 
occur significantly after the EOB light-ring [23] . 



PN accurate EOB waveform which already includes the 
leading NQC phase correction) and that it could simply 
be set to b = 0. We shall also do so here. 



IV. SELECTING BEST-FIT EOB FLEXIBILITY 
PARAMETERS 

As recalled in the Introduction, and in the previous sec- 
tion, the only EOB flexibility parameters whose best-fit 
values are still quite indeterminate are a 5 , v po i e and orr. 
In this section we shall show how to remedy this situation 
by combining information coming from various NR data, 
namely, on the one hand, from published Caltech-Cornell 
data, and, on the other hand, from recently computed 
Jena data (reported here). 



A. Using Caltech-Cornell published data to 
determine t> po i c and orr as functions of a 5 

To start with, let us recall that Ref. [21] had fixed 
Srr = and had then showed that imposing one con- 
straint relating the EOB waveform and Caltech-Cornell 
inspiral data, namely p^ d (a.5 7 w po ic) = 1, (see Eq. (35) 
in [21]), implied a rather precise functional relationship 
between w po i e and as (see Fig. 3 there). More recently, 
Ref. [22] compared the same type of EOB waveform with 
NR waveforms, computed with the CCATIE code of the 
Albert Einstein Institute, and suggested that it might 
be useful to flex the EOB waveform by introducing a 
nonzero value of <zrr, i.e. a non quasi-circular correct- 
ing factor / N< ^ c , Eq. (9), in the radiation reaction. Here 
we shall combine these two strategies by starting from 
an EOB waveform depending on the three a priori inde- 
pendent parameters (a 5 , 7j po io, orr) and by imposing two 
independent constraints relating the EOB waveform to 
published Caltech-Cornell data. These constraints have 
the form 



PoJ4 4 {05, flRR, fpolc) — 1, 

St' . , 
poj 4 (a5,aRR,i> P oio) = 1, 



(16) 
(17) 



where 



*T4EOB (*NR + St 



Pu m (a5,a RR ,v po i c ) = — 

<J m 

(18) 

Here A"™ 0t4eob (*nr + ) i s the value at the time 
*nr + ^*w m 01 the curvature waveform (tA 22 ) phase dif- 
ference between T4 and EOB when the two waveforms 
are matched at the gravitational wave frequency uj m (in 
the sense of [28]). The (5 m 's appearing in Eq. (18) are 
estimates of the value of the phase difference between 
TaylorT4 and Caltech-Cornell numerical relativity data 
at various times differing from the matching time i^g by 
6t LUm , as measured by us on the left panel of Fig. 19 
of [28]. Following the procedure outlined in Sec. IV 
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FIG. 3: Functional relationships linking « po i e and Srr to a 5 
obtained by imposing the two constraints (16)-(17) based on 
published Caltech-Cornell inspiral waveform data. 

of [21], we use the matching frequency iv 4 = uj m = 0.1. 
Then, we consider two of the measured values which 
have been used to produce the empty circles appearing 
in Fig. 5 of [21]. These two values are 

S 4 = 0.055 corresponding to 6t^ 4 = -1809M, (19) 
5' 4 = 0.04 corresponding to 5t' U4 = -529M. (20) 

The data point (5t UJ4 , 64) corresponds to the leftmost 
empty-circle on the top panel of Fig. 5 of [21], while the 
point {5t' bJ4 ,5' 4 ) corresponds to the next to next empty 
circle on the right of (<5^ 4 ,<$4). The former data point 
was used in Ref. [21] as the "main backward" W4 data. 
Note that the new data point that we use here is also 
"backward" (with respect to uo m = u 4 = 0.1), though 
it is less "backward" by about a factor three. We use 
these two points here because we think they represent the 
best "lever arms" to exploit the approximate 14 numerical 
data represented in Fig. 5 of [21]. In particular, we do not 
use any "forward" data point because the accuracy with 
which we could measure them is more uncertain. Let us 
emphasize that, as a consequence of this choice, our de- 
termination of the functional relationships v po i c (a^) and 
cirrus) exhibited below only relics on Caltech-Cornell 
data up to gravitational wave frequency Mu> < 0.1. 

The two constraints (16)-(17) were solved by numerical 
Newton- Raphson iteration in v po i starting from a grid of 



As a measure of the accuracy of the approximate data points 
quoted in Eqs. (19)-(20) above, let us mention that we have, 
since, directly determined from the Caltech-Cornell numerical 
data provided to us the values S4 = 0.05497 corresponding to 
8t Wi = -1809M, and 8' 4 = 0.03957 corresponding to St'^ 4 = 
-529M (with a numerical relativity time t„ 4 = 3782.1489M 
corresponding to u> 4 = 0.1). 



11 

TABLE II: Explicit values of the EOB effective parameters 
<2rr and t)poie for a certain sample of as. These values cor- 

.st' _. 
respond to imposing the two constraints p W4 ~ 1 ± 10 ~ 
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values of (orr, a§). The iteration was stopped when the 
constraints were satisfied to better than the 10~ 4 level. 
The result of this procedure consists of two separate func- 
tional relations linking, on the one hand, -y p oio to a 5 and, 
on the other hand, orr to a§. These two functional rela- 
tions are plotted in Fig. 3. The upper panel of the figure 
is a modified version of the v po \ c (a 5 ) functional relation- 
ship represented in the upper panel of Fig. 3 of [21]. The 
latter curve was drawn by fixing 5rr to zero and by im- 
posing only the first constraint, Pw 4 " 4 (a 5 , 0, v po ic) = L 
By constrast, the curve v po i e (ac,) in the upper panel of 
Fig. 3 was obtained by simultaneously tuning w po i and 
orr so as to satisfy the two constraints (16)-(17). The 
numerical data behind the plots of Fig. 3 are also given 
in explicit numerical form in Table II. [The many digits 
quoted there are only given for comparison purposes.] 

In the upper panel of Fig. 4 we exhibit, for the partic- 
ular value 125 = 25 (and, correspondingly, 5rr = 27.9197 
and Wp O i = 0.51563) the near-perfect agreement between 
the two W4 — matched phase differences 4>T4 — ^eob and 
(pT4 — 0nr- [Our choice of the particular value a§ = 25 
will be motivated in the next subsection]. For com- 
pleteness, we have also included in the upper panel (see 
dash and dash-dot curves) the analogous phase differ- 
ences matched at the matching frequencies L02 = 0.05 
and uj 3 = 0.063 instead of ^4 = 0.1. The visual agree- 
ment between these three phase-difference curves and the 
corresponding ones displayed in the left panel of Fig. 19 
in Ref. [28] is striking. [As in Fig. 19 of [28], we use here 
TaylorT4 3.5/2.5; see Appendix B for its precise defi- 
nition]. The lower panel of Fig. 4 plots the ^-matched 
phase difference (/)eob-0nr = [0T4-0nr]-[<At4-<Aeob], 
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Mw = 0.1 (in the sense of imposing the two constraints 
Eqs. (16)-(17)) we have reduced the number of indepen- 
dent unknown EOB flexibility parameters to only one, 
namely the "4 PN" EOB parameter 05. The basic phys- 
ical reason behind the difficulty of determining 05 by 
means of inspiral data only (especially when relying, as 
we did above) on data below GW frequency 0.1, is the 
fact that 0,5 starts significantly affecting the EOB dynam- 
ics (and waveform) only during the late inspiral, when the 
dynamics becomes strongly nonadiabatic. Our next step 
will be to constrain a$ by best fitting the EOB wave- 
form to numerical data covering more of the late-inspiral 
dynamics. 



B. Using numerical data covering late-inspiral and 
plunge to constrain the "4 PN" EOB flexibility 
parameter as 



FIG. 4: Top panel: near-perfect agreement between T4-EOB 
and T4-NR phase differences when as = 25, orr = 27.9197 
and fpoio = 0.51563. Here NR refers to the published results 
of the Caltech- Cornell inspiral simulation. The corresponding 
EOB-NR phase difference (bottom panel) is of the order of 
10 -3 radians over the 30 GW cycles of the Caltech-Cornell 
inspiral simulation. 



i.e., the difference between the two solid curves (red on- 
line and black) in the upper panel. 

Note that this phase difference varies between —0.0014 
and +0.0008 radians over the time span (of ~ 2464M up 
to Mui = 0.1) which was used in our EOB-CC fitting 
procedure. A study of the continuation of the curve ex- 
hibited in the bottom panel of Fig. 4 then shows that, af- 
ter Muj = 0.1 and up to a final frequency Mcu = 0.1565, 
this phase difference varies between +0.0013 and —0.0185 
radians over a time span of 96M. Note that a recent re- 
port [60] has indicated that a refined estimate of the total 
phasing error in the Caltech-Cornell simulation was of the 
order of 0.01 radians over the entire span of the simula- 
tion. Therefore the accuracy of our EOB-fit is consistent 
with such an error estimate. 

Summarizing so far: by best fitting the three- 
parameter flexed EOB waveform *f 2 OB (a5, v po ie, orr; t) 
to published 15 Caltech-Cornell inspiral data before 



5 Since we had had recently access to the actual Caltech-Cornell 
data we could and did check the reliability of the results ob- 
tained from the published data. In particular, when computing 
the phase difference <^bob — ^cc"* 1 we essentially recovered the 
results quoted in the text. For instance, we find that the ac- 
tual phase difference varies between: —0.002 radians at Caltech- 
Cornell time 600M and —0.01766 radians at the end of the sim- 
ulation (Mu> = 0.1565), passing through zero at Mw = 0.1. The 
number of GW cycles between tec = 600Af and tec = 3782 
(Mw = 0.1) is 22.10, while the number of GW cycles in the final 



In this subsection we shall fulfill, at least in first ap- 
proximation, the aim of the EOB-NR comparisons ini- 
tiated in Refs. [16, 20-22]; i.e., to determine an essen- 
tially unique set of "best-fit" EOB flexibility parameters 
(Upoie, orr ,, 05). In view of the results of the previous sub- 
section, we now need to best-fit the one-parameter flexed 
EOB waveform 

tf E 2 OB (a 5 ; t) = ^ B [v polc (a 5 ),a RR (a5),a 5 ; t], (21) 

where w po ie(o5) and cjrr^s) are the functional relation- 
ships illustrated in Fig. 3 above, to a numerical waveform 
smoothly connecting, without interruption, the nonadia- 
batic late-inspiral to the early-inspiral and to the subse- 
quent plunge. Here we shall make use of recently com- 
puted numerical data (see Sec. II) that cover (for the 
equal mass case) about 20 GW cycles of inspiral and 
plunge up to merger. As we shall see, for the purpose of 
determining ag, we will mainly use the signal only up to 
the plunge. 

As quantitative measure of the EOB-NR agreement 
we shall consider here the following norm of the 
«5-dependent EOB-NR phase difference (using the EOB 
metric waveform, Eq. (21) above) 

||A</C OBNR (a 5 ; t u t 2 ; t L ,t R ) = 

sup te[tLjtR] |0f 2 OB (a 5 ; t + Ti 2 ) + ai2 - <^ 2 R (g| . (22) 

Here [ti, , t R ] denotes the time interval on which one com- 
putes the Loo norm of the phase difference. In addition, 
(ti,t 2 ) denote two "pinching" times which are used to 
determine some time and phase shifts, Ti 2 = T(fi,£ 2 ) 
and a\2 = a(<i,£ 2 ), needed to compare the EOB and 
NR phase functions (which use different time scales and 
phase references). 



part of the simulation (after frequency 0.1) is 2.16. 
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FIG. 5: Loo norm of the EOB-NR late-inspiral ([£l, tn]) phase 
difference, as a function of as for v = 0.25 (1:1 mass ratio) 
and v ~ 0.2222 (2:1 mass ratio). NR refers to results of Jena 
coalescence simulations reported here. 



TABLE III: Pinching NR times and corresponding NR gravi- 
tational wave frequencies used to perform the EOB-NR com- 
parison of Fig. 6. 



V 


tx 


t 2 


22 


22 


0.25 


1764.9 


1940.1 


0.0998 


0.4716 


0.2222 


893.9 


1071.9 


0.1005 


0.4542 


0.16 


1297.6 


1476.3 


0.1051 


0.4189 



Let us recall the "two-pinching-times" procedure, in- 
troduced in [22], for determining the time and phase 
shifts t and a. First, the two waveforms being com- 
plex numbers, we decompose them in amplitude and 
phase: ^22 = Ay ex P(~ ^4> X ) where the label X can 
be either "EOB" or "NR" . The corresponding instan- 
taneous (metric) GW frequencies are then defined as 
uj x (t) = d(f) x /dt. We start by fixing two "pinching" 
times (£1,^2) on the NR time scale t. We then define the 
time-shift r by solving the equation </> NR (t 2 ) — </> NR '(£1) = 
EOB (r + i a ) - EOB (r + ti). Then, we define the phase 
shift a such that NR (ti) = EOB (ii + r) + a. In the 
limiting case where the corresponding GW frequencies 
u>x = cj nr (£i) and LO2 = a-> NR (i2) are nearly coincident, 
uji ps uj m k u>2, this procedure coincides with the one in- 
troduced in Ref. [28] and based on the choice of a single 
matching frequency uj m . 

We shall first consider the equal-mass case, v = 1/4. 
For this case we choose the following NR pinching times: 
ii = 1764.9 and t 2 = 1940.1 (corresponding to NR 
gravitational wave frequencies uif 2 = 0.0998 and w 22 = 



0.4717). These times bracket the merger time. This is 
done to optimize the EOB-NR agreement over the phys- 
ically most crucial (and possibly numerically most accu- 
rate) part of the waveform, i.e. the late-inspiral, plunge, 
merger and ringdown. Concerning the choice of the in- 
terval [^Lj ^r] used to compute the L m norm, we selected 
it with the following criteria in mind: as 05 is most im- 
portant during late-inspiral and plunge, but is somewhat 
uncorrelated to the way EOB approximates the plunge- 
ringdown matching, we chose [£l, £r] to cover the crucial 
stage of the late inspiral. More precisely, we have fixed 
ijt such that the NR gravitational wave phase is approx- 
imately 7.6 radians (i.e. 1.21 GW cycles) smaller than 
the phase when the EOB waveform modulus reaches its 
maximum (which is close to merger time in view of the 
discussion of Sec. III). Then, £l was chosen such that 
</> NR (£l)22 = 4>22"(t~R) ~ 58.5. Their numerical values are 
t L = 1198.8M and t R = 1899.6M, while the correspond- 
ing NR gravitational wave frequencies are W 22 = 0.05952 
and wj| 2 = 0.1898 respectively. Using these specified val- 
ues we have computed the L m norm of the EOB-Jena 
phase difference, Eq. (22), as a function of a^. The result 
is plotted, as a solid line (1:1 mass ratio), in Fig. 5. This 
figure shows that the limited range of values 20 < a 5 < 30 
is preferred in that it yields a minimum of the largest 
EOB-NR phase difference ||A<tj|| eobnr over the [£ l ,*r] 
interval specified above. This minimum phase difference 
is on the order of 0.01 radians. We note, in passing, that 
this late-inspiral interval partially overlaps (frequency- 
wise) with the range of the published Caltech-Cornell 
data as we used it above (i.e., focusing on frequencies 
Mu> < 0.1), but crucially extends to frequencies reach- 
ing roughly as high as the EOB adiabatic LSO frequency 
Olsc? = 0.2114). Though Fig. 5 is qualitatively simi- 
lar to the Loo norm of the EOB/Caltech-Cornell phase 
difference displayed in Fig. 4 of Ref. [21], it is impor- 
tant to remark that in the latter figure the L^ norm 
varied by only about a factor 2 over the entire as range, 
< CL5 < 100. By contrast, in the current Fig. 5 the 
Loo norm varies by about a factor 2 in the much smaller 
interval 15 < 05 < 35 and then increases by almost a 
factor 10 over the entire 05 range, 5 < 05 < 75. We 
can now use the "uncertainty level" ±0.026 radians in 
(f>2 2 R = 022 na (determined in Sec. II C above by compar- 
ing it with Caltech-Cornell data), as indicated by the 
horizontal line in the figure, to determine a correspond- 
ing interval of "best-fit" values of as . Though this uncer- 
tainty level is admittedly rather uncertain at this stage, 
it suggests that the "real" 16 value of 05 probably lies in 
the interval 12 < 05 < 40. To firm up our conclusion, 



Note that though G15 is, to start with, a theoretically well defined 
quantity within the EOB framework, its "experimental measure- 
ment" obtained by comparing specifically resummed versions of 
the EOB waveforms with numerical data partially transforms it 
into an "effective parameter" describing a complicated nonper- 
turbative process. 
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we have also considered numerical data concerning the 
2:1 mass ratio case. In that case we considered again 
the Loo norm, Eq. (22), and we made similar choices 
both for the pinching times and for the extremities of 
the Loo interval. In particular, £r was chosen to sit 7.6 
radians before the maximum modulus while we kept the 
left-right phase difference to the same value as above, 
namely ^ 2 R (*l) = ^ 2 R (i R ) - 58.5. The resulting LooK) 
function is plotted as a dashed line in Fig. 5. Though 
the minimum of this curve is much more shallow than 
before, the important fact is that the 1:1 preferred as 
range is consistent with the 2:1 result. Let us ob- 
serve (without wishing to attribute any deep significance 
to this fact) that the preferred range for 05 happens to 
be close to the "special" (15 value for which the "EOB- 
horizon" decreases, when v increases up to 1/4, down to 
a vanishing EOB radial coordinate. Indeed the P\ Pade 
approximant that we use here to define the as-flexed EOB 
radial potential A{u) = P\ L4 Taylor (u)] has the structure 
A(a$; u) = (1 — rHU)/Di{u) where D^u) is a 4th-order 
polynomial in v (see Eq. (3.10d) of Ref. [17]), and where 



768 - (3584 - 123tt> - 2Aa b v 
1536 - (3776 - 123tt 2 )u ' 



(23) 



Here th is the radial location of the "EOB horizon" , in 
the sense that A(u) vanishes for r = 1/u = rjj (at least 
when Th is positive). For any given positive 05, ru is 
a decreasing function of v. If we require that rjj stays 
positive for all values of v £ [0, 1/4], we find that 05 must 
be smaller than the "special" value 



123tt 2 - 512 
24 



29.2484. 



(24) 



Note, however, that there is nothing a priori wrong with 
higher values of 05. In that case the radial function A(r), 
considered versus r, has anyway a third-order zero at 
r = 0. 

Summarizing: by combining the comparison of the 
EOB waveform with, on the one hand, published Caltech- 
Cornell inspiral data and, on the other hand, our coa- 
lescence data, we have been able to select a preferred 
small region of the EOB flexibility parameters. This 
region is made of (approximately) correlated triplets 
(05, v po \e(a5), orr^)), and is located between the sec- 
ond and the 8th lines of Table II. 



V. DETAILED EOB-NR WAVEFORM 
COMPARISONS FOR a 5 = 25 

To confirm the validity of the conclusions reached in 
the previous section, we shall now study in detail the 
performance of the center of the above selected interval, 
namely 0,5 = 25 together with the corresponding values 
of orr and v po i c listed in Table II. 

In this section we shall consider numerical waveforms 
for three different values of v, namely v = 0.25, v = 2/9 = 



0.2222 and v — 0.16 (corresponding respectively to the 
mass ratios 1:1, 2:1 and 4:1) extracted from the simula- 
tions of Table I. Note that our best-fit procedure out- 
lined above essentially relied only on the 1:1 mass-ratio 
case so that the other cases that we consider here will 
test the ability of the EOB formalism to capture the NR 
waveforms. The EOB flexibility parameters used for the 
various mass ratios are the ones listed in the fifth row of 
Table II. In view of the proximity of the "best-fit" i> po i c 
value v£fe(v = 0.25) = 0.5156 to the "best-fit" v po i c 
found (following the strategy of [21]), in the test mass 
limit, v^l(f = 0) = 0.52655 (for the P 4 4 4 PN-accurate 
flux), we made no attempt at interpolating v po ie(^) be- 
tween the two values of v. 

To compare EOB and NR waveforms we follow the pro- 
cedure indicated above. This procedure involves choos- 
ing two "pinching" times t\ and t 2 (which should not be 
confused with the times t\> and £r which will play 
no role in this section). We summarize in Table III the 
"pinching" times we use, together with the corresponding 
frequencies. Note that in all cases the lowest pinching fre- 
quency is around 0.1 while the highest one (reached after 
the merger) is roughly 10% lower than the main ringdown 
frequency. 

The results of the detailed EOB-NR comparison arc 
presented in Fig. 6 and 7. For completeness, we have 
used the full numerical waveforms including the burst of 
junk radiation it contains at the beginning. 

The two upper panels of Fig. 6 refer to the equal-mass 
case [y = 0.25). On the left, we plot the "pinched" EOB- 
NR phase difference (in radians) over the full simulation 
time (see inset). Note that the total simulation covers 
about - 146 radians of GW phase; i.e., 23.24 GW cy- 
cles (starting from the beginning of the inspiral, when 
£nr ~ 110M, to the middle of the ringdown, up to 
i N R = 1980M). We see that the EOB-NR phase dis- 
agreement stays quite small during most of the inspiral. 
More precisely A0 EOBNR stays in the range [-0.04, 0.01] 
all over the time interval 1200 < iNR S 1930. This cor- 
responds to a "two-sided" (in the sense of footnote 12 of 
Ref. [22]) EOB-NR phase difference smaller than ±0.025 
radians, or ±0.004 GW cycles over 730M. As in pre- 
vious analysis, the jump in the phase difference around 
£nr ~ 1930 is connected to the rather coarse way in 
which the EOB formalism represents the merger. Still, 
the accumulated phase difference over the transition be- 
tween plunge and ringdown is only of the order of 0.15 
radians; i.e., 0.02 GW cycles. Note that over the full 
simulation time (see inset in top-left panel) there is an 
accumulated phase difference of about -0.2 radians. In 
view of the discussion on the accuracy of the numeri- 
cal simulations in Sec. II, it is quite possible that this 
difference is mainly due to effects related to the use of 
finite extraction radii. Similarly, part of the phase dis- 
agreement around the merger might come from numerical 
inaccuracies. The upper right panel of the figure com- 
pares the real part of the two metric waveforms. The 
visual agreement between the two is striking, apart from 
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FIG. 6: Comparison between NR and EOB waveforms for v = 0.25 (top), v = 0.2222 (bottom). The left panels depict the 
EOB-NR phase difference; the right panels show the real part of the metric waveforms. Here, NR refers to the full results of 
Jena coalescence simulations, from early-inspiral to ringdown, by contrast to the Loo norm of Fig. 5 which concerned a late 
inspiral stage [tL,tit]. This interval is indicated on the figures. The pinching times (£i,ta) of Table III are also shown. The 
dash-dot and the dash vertical lines at the extreme right of the figures mark the location of the EOB adiabatic LSO and the 
"EOB light-ring" respectively. 



the amplitude disagreement (~ 20%, see below) localized 
around the merger. In view of the discussion in Sec. II, 
part of this difference might also have a numerical origin. 

The bottom panels of Fig. 6 refer to the 2:1 mass ra- 
tio case {y — 2/9 = 0.2222). Here the phase agreement 
(left panel) is even better than before. Over the nearly 
full time interval 143 < t NR < 1100 the EOB-NR (two- 
sided) phase difference is smaller than ±0.05 radians; 
i.e., ±0.008 GW cycles. The corresponding middle-right 
panel compares the real part of the two metric waveforms. 
Again, the agreement is striking apart from a ~ 20% am- 
plitude disagreement localized around the merger (see 
below) . 

Finally, Fig. 7 deals with the 4:1 mass ratio case (v = 



0.16). Here the agreement is still quite good, though it is 
noticeably less good than in previous cases. Consistently 
with the discussion of numerical accuracy in Sec. II, this 
less compelling accordance is likely to have its origin in 
numerical discretization errors. A clarification of this 
issue would need higher-accuracy simulations. 

Figure 8 completes the comparison between EOB and 
NR waveforms, for the equal-mass ratio case, by simulta- 
neously displaying, versus time: (i) the two GW frequen- 
cies 17 , (ii) twice the EOB orbital frequency O, and (iii) 



For clarity we add in several figures a subscript 22 to the gravi- 
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FIG. 7: Comparison between NR and EOB waveforms for v = 0.16. The left panel depicts the EOB-NR phase difference; the 
right panels show the real part of the metric waveforms. The pinching times (t\,to>) of Table III are also shown. The dash-dot 
and the dash vertical lines at the extreme right of the figures mark the location of the EOB adiabatic LSO and the "EOB 
light-ring" respectively. 



the two moduli. The leftmost (dashed) vertical line indi- 
cates the location of the EOB adiabatic LSO, while the 
rightmost one refers to the "EOB-light-ring" . Though 
this figure exhibits the approximate nature of the EOB 
matching procedure (notably visible in the small differ- 
ences in the GW frequencies), it also illustrates how the 
apparently coarse EOB-matching procedure is able to ef- 
fectively reproduce, with high accuracy, the overall time 
variation of the GW frequency through the merger onto 
the ringdown. We have obtained similarly good agree- 
ments for the other mass ratios. 

We conclude this section by showing in Fig. 9 the 
fractional amplitude differences, for the three mass ra- 
tios considered here, between EOB and NR waveforms. 
The solid line in the figure plots the quantity AA/A = 
(Aeob — -AnrV-^nr versus NR time for v = 0.25. It is 
quite possible that the approximately linear trend vis- 
ible on this (solid) line is due to effects related to the 
finite extraction radius; the decrease in amplitude dis- 
agreement as we go to later inspiral times is consistent 
with the decrease in amplitude uncertainty (as discussed 
in Sec. II) as the amplitude rises. If this is the case, the 
minimum value, before the merger, might be indicative of 
the actual EOB-NR amplitude agreement. For v = 0.25 
this minimum is min[AA/A] rj +5 x 10~ 3 . The jump 
in A A I A during merger is of the order of 20%. Though 
part of this jump might have a numerical origin, we think 
that most of it comes from the EOB approximate match- 



tational wave frequency or phase as a reminder of the fact that 
we compare quadrupolar t = m = 2 waveforms. 



ing procedure around merger. Let us recall, in this re- 
spect, that in Ref. [22] AA/A, for v — 0.25, was of order 
±1% during inspiral and rose to a maximum of +18% at 
merger. The leftmost curve (dashed line) on Fig. 9 refers 
to the v = 0.2222 case, while the middle curve (dash-dot 
line) refers to the v — 0.16 case. For the same reasons 
as above it is likely that the approximate linear trends 
(which are smaller by a factor ~ 2 than before) are of 
numerical origin. The minimum values before merger of 
AA/A are min[AA/A] w +7 x 10~ 3 (for v = and 
min[AA/A] » +5 x 10~ 3 (for v = 0.16). Note that the 
jumps in [AA/A] around merger are quite similar to the 
v = 0.25 case, namely about ~ +20%. 



VI. CONCLUSIONS 

We have compared the "flexed" [17, 18] resummed 
3 +2 PN-accurate [21] Effective-One-Body (EOB) wave- 
form to two, independent, numerical relativity (NR) data 
on inspiralling and/or coalescing binary black hole sys- 
tems: on the one hand, published Caltech-Cornell inspi- 
ral data [28] (mainly used by us only up to Muj < 0.1) 
and, on the other hand, newly computed coalescence data 
using the BAM code [30, 31]. 

Wc effected this EOB-NR comparison with a strategy 
allowing us to locate a "best-fit spot" in the space of the 
three main EOB flexibility parameters (05, v po i e , orr). 
This strategy is multi-pronged: 

• We selected two measurements of published 
Caltech-Cornell equal-mass inspiral data concern- 
ing the TaylorT4-NR phase differences at two dif- 
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FIG. 8: Comparison between EOB and NR instantaneous 
gravitational wave frequencies (and moduli) for the equal 
mass case, v = 0.25. Here, as in Fig 9, NR refers to the 
Jena coalescence simulation. The dash-dot vertical line indi- 
cates the EOB adiabatic LSO, while the dash one the "EOB 
light-ring". The pinching frequencies (u)\ 2 \ui2 2 ) of Table III 
are also indicated. 



ferent times, approximately spanning the GW fre- 
quency interval 0.04 < ui < 0.1. 

• We imposed two constraints requiring that these 
NR phase differences be equal to two correspond- 
ing analytical TaylorT4-EOB phase differences, see 
Eqs. (16)-(17). This gave us two equations for the 
three main flexibility parameters (05, v po i c , orr). 
By numerically solving these two equations we de- 
termined two functional relationships linking, sep- 
arately, fpoic to 05 and 5rr to 05. See Fig. 3 and 
Table II. 

• Having in hands these "Caltcch-Cornell-prefcrred" 
functional relationships w po i e (a5) and cirrus), we 
selected from our newly computed coalescence sim- 
ulation (again for the equal-mass case ls ) a time 
interval [£l,£r] corresponding to the following GW 
frequency interval 0.060 < uj < 0.19. On this time 
interval we compared the Jena numerically simu- 
lated phase evolution to the as-dependent analyt- 
ical EOB one, and we computed the L m norm of 
their difference, i.e. (see Eq. (22) for more details) 



l|A0|| 



EOBNR 



(05; tL, £r) 



su Pte[t L ,tR] \4>22 B (a5; t) 



(25) 



FIG. 9: Fractional EOB-NR differences in the gravtitational 
wave (metric) amplitudes, A = ^22!, versus NR time for 
different mass ratios. 



• We plotted (as a solid line) in Fig. 5 ||A0||EOBNR 
as a function of (25 . We found that this norm 
has a rather well localized minimum around the re- 
gion 20 < 05 < 30. To transform this fact into an 
actual "error-bar" on the value of 05 we would need 
to have in hands a precise measure of the level of 
the errors present in the (Jena) numerical data over 
the time interval [finfe] on which the norm is 
computed. At this stage we do not have at our 
disposal a reliable measure of this error. However, 
in Sec. II C we have given what we think is our 
current best estimate of this error level by directly 
comparing, on the crucial time interval [£l, £r], the 
Jena phase data to the actual Caltech-Cornell data. 
This current best estimate is ±0.026 radians and, 
according to Fig. 5, would correspond to the fol- 
lowing confidence interval for 05: 12 < 05 < 40. 
More work is needed to nail down in a more precise 
way the error level in the Jena phase (see in par- 
ticular our discussion above on the internal error 
estimate based on comparing various radius extrap- 
olations methods). In addition, for any value of 05 
in such an allowed confidence level, we would con- 
clude that the corresponding triplets of correlated 
values 05, Wpoie(a5) and <2RR(a5) obtained from Ta- 
ble II, determine preferred best-fit values of the 
EOB flexibility parameters 19 ((15, v po i c , 5rr). In 



When best-fitting the EOB flexibility parameters, we use v = 
0.25 data because these are more sensitive to 015. 



Note that the functional relationships ■u P oie(i5) and ajm(as) dis- 
cussed above have no invariant physical meaning and are just 
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other words, our current preferred values of the 
EOB parameters (a^, v po i e , (Xrr) lie between the 
second an the 8th lines of Table II. 

• The present implementation of this strategy is, 
however, certainly somewhat affected by numeri- 
cal noise. A possible indication of this fact is that 
the computation of a similarly selected norm 
pertaining to the 2:1 mass ratio simulation gives 
results (plotted as a dashed line in Fig. 5) which, 
though they are fully consistent with the 1:1 mass 
ratio case, exhibit a more shallow minimum versus 
a 5 . For the 2:1 diagnostic to select an inter- 
val of preferred values of a 5 we would need a reli- 
able estimate of the numerical error level in the 2:1 
phase data. However, at this stage we do not have 
such an estimate. The rough error level quoted 
in Sec. II D is just a very conservative upper limit 
which, moreover, does not concern the specific time 
interval [£l,£r] wc arc interested in. Let us em- 
phasize that, anyway, even if as is allowed to vary 
in the full interval 5 < a 5 < 75 that we explored, 
the maximum EOB-NR phase disagreement (on the 
considered late-inspiral interval, which corresponds 
to about 58.5 radians before crossing the last stable 
orbit) is below 0.1 radians, i.e. 0.015 GW cycles. 

• We think that it would be necessary to devote 
a special effort toward having very high-accuracy 
numerical simulations covering the crucial late- 
inspiral, corresponding to the frequency range 
0-1 ^ w ~ 0-2; for several mass ratios. Pending the 
availability of such simulations, wc provisionally 
conclude that our current "best-bet" choice of EOB 
flexibility parameters is at the center of the above- 
selected interval: i.e., it is given by 0,5 ~ 25 together 
with the correlated values of u po ic and orr listed in 
Table II. In Sec. V we presented evidence that these 
values of (05, v po i e , 5rr) lead to an excellent agree- 
ment between EOB and NR for several mass ratios 
and for the entire time-interval covering inspiral, 
late-inspiral, plunge, merger and ringdown. In par- 
ticular, we found that the dephasing between EOB 
and our new coalescence data are smaller than: (i) 
±4 x 10~ 3 GW cycles over 730M (11 cycles), in the 
equal mass case, and (ii) ±8 x 10~ 3 GW cycles over 
about 900M (17 cycles) in the 2:1 mass-ratio case. 
In addition, we recall that the phase difference be- 
tween our current "best-bet" EOB and both pub- 
lished and actual Caltech-Cornell data stays within 
0.018 radians over the entire span of the simula- 
tion. Such a phase inaccuracy is comparable with 
the current, updated estimate of the numerical er- 
rors of the waveforms of Ref. [28], namely 0.01 ra- 



intermediate tools in converging on the looked-for best-fit point 
in the three dimensional EOB flexibility parameter space. 



dians [60]. 

• As a contrast to the EOB performance, we also 
study in Appendix B the performance of the Tay- 
lorT4 approximant. Our analysis shows that the 
apparently good performance of TaylorT4 during 
the inspiral is due to a lucky compensation between 
two effects going in opposite directions: (i) the bad 
convergence of the adiabatic PN expansion and (ii) 
the fact that the T4 approximant does not take into 
account nonadiabatic effects. This compensation 
causes an "enhancement" in the domain of validity 
of T4. However, we show that this enhancement 
holds only for a limited range of values of the mass 
ratio. This is consistent with the finding of [65] 
that the enhanced validity of T4 is fragile and is 
undone by spin effects. 

In conclusion wc think that the results presented here 
corroborate the aptness of the EOB formalism to pro- 
vide accurate representations of general relativistic wave- 
forms. We suggest that the specific 3 +2 PN-accurate re- 
summcd EOB waveform (with the current "best-bet" 
values of the flexibility parameters determined here) be 
used in constructing banks of waveform templates for cur- 
rently operating gravitational wave detectors. 
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APPENDIX A: COMPUTING METRIC 
WAVEFORMS FROM CURVATURE 
WAVEFORMS 

This first Appendix is devoted to the discussion of an 
appropriate way of choosing the integration constants 
that enter the metric waveform h(t) when deriving it by 
double time-integration from a given (numerical) curva- 
ture waveform ip4(t). 

Our conventions are as follows: for reasons of continu- 
ity with the recent papers [18, 22, 66] we use the normal- 
ization factor N e = ^/(i + 2)(£ +!)£(£ - 1) in the metric 
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waveform to get the so-called Zerilli-Moncrief normal- 
ized waveform that we shall denote *^ o) (for even and 
odd-parity) as used in Ref. [67] . The metric waveform is 
expanded in spin-weighted spherical harmonics of spin- 
weight s = — 2 as 

oo i 

h+-ih x =J2Yl h im -2Y tm (8, 4>) (Al) 

1=2 m=-£ 

where the link between the multipolar metric waveform 
hl m (as used for instance in [68]) is 

tf" = ^(*S + i*S). (A2) 

The raw output of the numerical simulation used here is 
the Newman-Penrose scalar ip^.. This is decomposed in 
harmonics as 

oo I 

h+ - ih x =V^4 = E E 4i m -2Y tm (0, <f). (A3) 

1=1 m=-l 

The computation of the Zerilli metric multipoles from 
its curvature correspondant ■0| m requires a double time 
integration. Various ways of fixing the two integration 
constants entering this process have been discussed in 
the literature [15, 62, 69, 70]. We focus here on the 
I = m = 2 multipole of the Zerilli-Moncrief normalized 
metric waveform . 

We wish to emphasize that the choice of integration 
constants is particularly delicate when dealing with the 
metric waveform h(t), by contrast to dealing with the 
quantity h(t) which is most prominent in other appli- 
cations, such as the computation of recoil. For instance, 
Ref. [70] suggested to integrate backward in time starting 
with zero integration constants at t = +00. This proce- 
dure leads to a rather accurate h(t). However, we found 
that the resulting h(t) is not accurate enough for the pur- 
pose of high-accuracy waveform comparison discussed in 
this paper. This is exemplified in Fig. 10. This figure 
shows the metric waveform obtained by such a backward 
integration. The important point is that the modulus of 
the complex waveform exhibits quite visible unphysical 
oscillations at early times (during inspiral). 

By contrast, we found that the following procedure (re- 
lated to some of the suggestions of Ref. [62] ) gave reliably 
accurate results. We start by computing (e.g., separately 
for the real and imaginary parts, or directly for the com- 
plex quantity) the first and second forward time integrals 
(using e.g. Simpson's rule) of rip^, starting at t = with 
zero integration constants, i.e., we define 

ho(t)= [ dt'ri4 m (f), (A4) 
Jo 

h (t) = f dt'h Q (t'). (A5) 
Jo 
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FIG. 10: Computation of the Zerilli normalized metric wave- 
form 9^ fr° m r " i / , l 2 y i a two backward time integrations start- 
ing with zero integration constants at the final time. The data 
refer to the 1:1 mass ratio (y — 0.25) numerical simulation. 
Unphysical oscillations in the modulus are quite visible at 
early times. 

Then, we fit over the full simulation time interval (sepa- 
rately for the real and imaginary parts) the second inte- 
gral ho(t) to a linear function oft, say /ig in-fit (t) = at + (3 
where a and j3 are complex quantities. Finally, we define 
the metric waveform as 

h(t) = h (t) - h% n - Rt (t) = h (t) - {at + 13) . (A6) 

Note that this also defines the time-derivative of the met- 
ric waveform as 

h(t) = ho(t) -a. (A7) 

The result of this procedure is shown in Fig. 11. Here 
we applied the procedure explained above to the ripj 2 
waveform coming from the 1:1 mass-ratio simulation ex- 
tracted at r — 90. The top panel shows the real and 
imaginary parts of ho(t) (divided by the normalization 
factor N2) together with their best linear fits, i.e. the 
real and imaginary parts of at + /3. The bottom panel 
shows the final waveform h(t), i.e. the difference between 
ho(t) and the best linear fit at + (3. The important point 
is to notice that the modulus of h(t) (the blue line in the 
bottom panel) is monotonically increasing with t during 
inspiral without exhibiting any of the unphysical oscilla- 
tions that were present in the previous figure. 20 We show 



This is a good indication that the integration constants have been 
computed correctly and that the real and the imaginary parts of 
the waveform are dephased by 7r/2 with very good approxima- 
tion. 
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FIG. 11: Computation of the metric waveform from ripf 2 
via two time integrations starting at t — 0. The upper panel 
has zero integration constants and exhibits clear linear drifts 
at + [3. The bottom panel shows the result of subtracting 
the linear drift of the waveform obtained by fitting the upper 
panel over the entire time interval starting at t = 0. 
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FIG. 12: Same as Fig. 11 except that the integration and the 
linear fit have been done starting at time t ~ 150, i.e. at 
the beginning of the inspiral signal. Note the oscillations in 
the modulus of the bottom panel, and the fact that the linear 
drifts (visible in the upper panel) are much larger than in 
Fig. 11. 



on the same plot also the real and imaginary parts of the 
complex quantity (which correspond to the h+ and 
h x polarizations of the wave after division by r and mul- 
tiplication by the spin-harmonic ^Y 22 ) as well as the 
gravitational wave frequency 0^22 obtained from the met- 
ric waveform an d the gravitational wave frequency 
obtained from the curvature waveform ip 22 . 

In addition, let us emphasize that for this procedure to 
work it is important to start the integration from the ab- 
solute beginning of the numerical simulation, by which we 
really mean t = 0, i.e. before any signal reaches the ob- 
server. One might have thought that it is better to start 
the integration after the junk radiation, at the beginning 
of the inspiral signal. This is not the case, as it is illus- 
trated in Fig. 12. This figure shows the worsened result 
we obtain when we use exactly the procedure explained 
above, but on the time interval t > 150, i.e. starting at 
the beginning of the inspiral signal instead of starting at 



t = 0. Note the oscillations in the modulus of \& 



(e) 



By 



contrast, even if we blow up the corresponding graph in 
Fig. 11 the oscillations are practically absent. Note also 
that the linear drifts are now much larger than before. 
This is part of the reason why the results are less good in 
this case. By contrast to the first case where, starting at 
t = meant starting with extremely small initial values 
of rip 22 , starting at the beginning of the inspiral means 
starting with much larger values of ri/j 22 : this effect en- 
larges the linear floors and therefore the errors on the 
determination of the linear floors to be subtracted. 



APPENDIX B: SHORTCOMINGS OF THE 
TAYLORT4 PN APPROXIMANT 



To contrast with the EOB-NR comparison done in the 
text, we consider in this appendix the comparison be- 
tween the so-called TaylorT4 post-Newtonian approxi- 
mant [23, 25, 28, 65, 71, 72]. and various NR data. 
This approximant is defined by two successive prescrip- 
tions: the first concerns the computation of a "T4 orbital 
phase" $T4(i) while the second concerns the definition 
of a "T4 metric waveform" . Here we shall focus only on 
the I = m — 2 quadrupolar waveform. The "T4 orbital 
phase" $T4(i) is defined by integrating the ordinary dif- 
ferential equations 



T4 



(It 

dx 



r 3/2 



64z/ 



5 Taylor 



dt 



-x a 



3.5 



(Bl) 
(B2) 



where is the 3.5 PN Taylor approximant, for any 

given value of v, to the Newton-normalized ratio (flux- 
function) /(derivative of energy hmctioii)—F(x)/E(x) 
where E'(x) — dE/dx. As in the text, we scale dimen- 
sionful quantities by the total "bare" mass M — m 1 +m 2 . 
This is for instance the case for the time variable t in the 
above equations. The explicit expression of a^ lm (x) 
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reads [23, 25] (for the nonspinning case) 
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This phasing evolution is completed by a quadrupolar 
waveform which is known (for any given value of v) at 
the 3 PN accuracy level [21, 62, 68, 73]. Following [28, 68] 
we define the 3 PN-accurate T4 waveform by dropping all 
the ln(x/x ) terms in Eq. (79) of Rcf. [68]. We display it 
explicitly here to clarify which waveforms we use in our 
T4 studies. The explicit expression of the I = m = 2 
Zcrilli-normalizcd metric waveform reads 
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where je = 0.57721 ... is Euler's constant. The taylorT4 
3.5/2.5 waveform (used in most of our comparisons) is 
obtained by dropping the terms oc x 3 [aln(a;) + b] on the 
r.h.s. of this equation. 

Thorough comparisons between the TaylorT4 3.5/2.5 
waveform (i.e., 3.5 PN accuracy for phase and 2.5 PN ac- 
curacy only for amplitude) and NR waveforms were per- 
formed, for the equal mass case, v = 0.25, in [28, 65, 71]. 
Ref. [28] concluded that this approximant yields an "as- 
tonishingly good" agreement with numerical data during 
the inspiral, i.e. a dephasing smaller than 0.05 radians 
over ~ 30 GW cycles before reaching the GW frequency 
Mu>22 =0.1. On the other hand, Ref. [65] showed that 
the inclusion of spins on the black holes had the effect of 
considerably worsening the agreement between T4 and 
NR data. Here we shall study the effect of varying the 
mass ratio (for nonspinnig black holes). We shall also 
go beyond the analyses of [28, 71] in discussing the be- 
haviour of T4 for GW frequencies above 0.1. Let us first 
compare the 21 v — 0.25 TaylorT4 3.5/2.5 quadrupo- 
lar waveform if/Ji with equal-mass NR waveforms com- 
puted by the Jena group. As discussed in Sec. II, the 
BAM code outputs the Newman-Penrose curvature scalar 
ip4(t,r 7 9 7 (p) at various extraction radii r. This angular- 
dependent curvature scalar is then: (i) decomposed on 
the basis of spin-weighted spherical harmonics and then 



Here, to facilitate the comparison with previous work, we use 
a T4 approximant with 2.5 PN accurate amplitude. Our main 
conclusions would be similar had we used the 3 PN accurate 
amplitude. 



(ii) integrated twice over time to yield the metric wave- 
form *22 R - The choice of integration constants in this 
integration procedure was done according to the proce- 
dure outlined in Appendix A. To compare the two wave- 
forms an d ^22 > as functions of their respective time 
variables, we choose a relative time-shift r and a relative 
phase-shift a by following the same two-pinching-time 
procedure used in the text. 

Figure 13 compares the gravitational wave frequency 
o>22 computed from the numerical data and plotted as a 
function of the NR time scale (solid line) with that of 
the TaylorT4 3.5/2.5 approximant plotted as a function 
of the shifted T4 time-scale (dash-dot line). The two 
waveforms have been "pinched" at the NR times t\ — 
1299.9 and ti = 1399.8, corresponding to NR frequencies 
cji = 0.062643 and uj 2 = 0.066292, respectively (which 
approximate the matching frequency io m = ^3 = 0.063 
of [28]). 

We see on this figure that there is a very good agree- 
ment between the two frequencies during the inspiral, up 
to, say, the NR time t NR = 1850, where w^ 2 R = 0.1301 
and L0J2 differ by about 2.2%. Then, soon after, ojJ 2 4 
starts deviating very strongly from and blows up to 
infinity at the finite time ibiowup = 1899.5 (indicated by 
the leftmost vertical dashed line in the figure). This blow- 
up time, considered on the NR time-scale, corresponds 
to a NR frequency w££ wup ~ 0.1889. Note that this fre- 
quency is smaller than the Effective-One-Body GW fre- 
quency at the adiabatic Last Stable Orbit (LSO), which 
is equal to 2£1ls'o = 0.2114 (corresponding to an EOB 
radial coordinate rf^o — 4.4729). The rightmost verti- 
cal dashed line in the figure indicates the "u> — LSO" , in 
the sense of Ref. [2], i.e. the time when the (NR) GW fre- 
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FIG. 13: Numerical relativity (Jena) and TaylorT4: compari- 
son between the instantaneous gravitational wave frequencies 
for the equal mass case (y = 0.25). 



quency w^ 2 R equals the adiabatic LSO frequency. Here we 
consider the case v = 0.25 and a$ — 25 and we compute 
the LSO frequency within the EOB approach. Therefore, 
in the equal mass case, the TaylorT4 approximant breaks 
down already during late inspiral, before the EOB LSO 
and before the plunge. 

The fact that the T4 approximant blows up at a fi- 
nite time is a simple mathematical consequence of the 
structure of the differential equation (B2), given that, 



a 



3 5 yl ° r ( a: ) i s f° un d to remain positive for every x > 0. In- 
deed, one can even easily analytically compute the blow- 
up time as being ibiowup = to + At where to is any given 
"reference" time on the T4 time scale (corresponding to a 
frequency parameter x(to) — Xq), and where At is given 
by the following convergent integral 



At = 



clx 



Xq OjyX 1*3 5 



(B6) 



where C v = 64^/5. 

After having compared the T4 approximant to NR 
data (in the equal-mass case) let us compare the T4 ap- 
proximant to the EOB one. As emphasized in Ref. [21], 
a convenient way of comparing two waveforms (which 
avoids the issue of finding suitable time shifts and phase 
shifts) consists in considering the following shift-invariant 
"phase-acceleration" function 
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Note that in the present paper we consider the frequen- 
cies of the metric waveforms (by contrast to the frequen- 
cies of the curvature waveforms considered in [21]). 

In the left-panel of Fig. 14 we compare the phase ac- 
celeration curves of T4 3.5/2.5 (dash-dot line) and EOB 



(solid line) for the equal mass case. The leftmost verti- 
cal line indicates the EOB frequency w 0.19 correspond- 
ing to the T4 blow-up time (computed by Eq. (B6)). 
The rightmost vertical line indicates the adiabatic EOB 
LSO frequency, 20Lg ) Q ~ 0.21 as above. We terminated 
the horizontal axis at UI22 = W22 = 0-3676 which cor- 
responds to the EOB time when the EOB orbital fre- 
quency reaches its maximum; i.e., the so-called "EOB 
light-ring" , which defines the "merger time" within the 
EOB approach. Note that this figure shows the metric 
waveform analogue of the EOB (curvature) phase accel- 
eration curve of Fig. 2 of Ref. [21] and extends it up to 
the merger time. As was already emphasized in [21], the 
figure shows that the T4 acceleration curve strongly di- 
verges away from the EOB one for frequencies LO22 ^ 0.1, 
i.e. during the late inspiral, before reaching the LSO. 

The right-panel of Fig. 14 illustrates the case where 
the mass ratio is 4 : 1, i.e. v = 0.16. The vertical 
dashed line indicates the adiabatic EOB LSO frequency 
2Qlso ~ 0.17. For this value of v, the blow-up frequency, 
computed as above, turns out to be larger than the EOB 
light-ring frequency LO22 — 0.3201. We see on this plot 
that, contrary to the equal mass case, the T4 acceleration 
curve starts to deviate significantly from the EOB one for 
frequencies W22 > 0.05. Note, however, that because the 
two curves cross again just before the LSO, we expect 
that the phase difference between T4 and EOB will re- 
main, on average, rather small up to the LSO. However 
later on the T4 phasing will drastically deviate from the 
EOB one. 

Finally, Fig. 15 considers the test-mass limit (y — > 0). 
Here we compare three acceleration curves: (i) the adi- 
abatic limit of the T4 acceleration curve, given simply 
by a u (u>) = a'^ OT (x) computed in the limit v — » 
and with x = (cj/2) 2 / 3 ; (ii) the EOB a w curve com- 
puted for v = 0.01 and (iii) the exact adiabatic limit 
of the test-mass acceleration curve, i.e. the Newton- 
normalized ratio F{x)/E'(x) (see e.g. Ref. [7]). Here 
the flux function F(x) is the one computed numerically 
in Ref. [74, 75]. The two vertical lines in the figure refer 
to the v = limit of the adiabatic LSO frequency (left- 
most line, 20 = 0.1361 ) and to the v = 0.01 EOB adi- 
abatic LSO frequency (rightmost line, 20lso = 0.1378). 
This figure illustrates two facts: first, the T4 approxi- 
mant starts strongly deviating from the exact result early 
on (say for 0J22 < 0.06, see inset); second, one needs 
to consider v < 0.01 to ensure that the usual adiabatic 
approximation is satisfactory up to frequencies close to 
the LSO one. This is consistent with the analytical 
estimate obtained in [2] according to which the devia- 
tions from adiabaticity become important when the fre- 
quency fractionally deviates from the LSO frequency by 
(Jw/wlso ~ ^ 2/5 - The presence of the 2/5 power means 
that we need v < 3 x 10~ 3 to be approximately adiabatic 
up to 90% of the LSO frequency. 

Summarizing, the main results of the present Appendix 
(and of complementary investigations of the different 
"speeds" with which the T4 and EOB waveforms "move" 
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- TaylorT4 3.5/2.5 

— EOB 





FIG. 14: Comparison between phase acceleration curves a u ol EOB (solid line) and TaylorT4 (dash-dot line) for v = 0.25 (left 
panel) and v = 0.16 (right panel). The rightmost dashed vertical line indicates the location of the EOB adiabatic LSO. The 
leftmost vertical line on the left panel indicates the EOB GW frequency corresponding the the instant when T4 blows up. 
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FIG. 15: This figure (done for v <C 1) illustrates, in particular, 
the fact that the domain of validity of TaylorT4 when v < 0.16 
reduces to the "normal" 3.5 PN one, namely u < 0.06. 

as v varies) are: (i) we predict that the T4 approximant 
will define an effective phasing template for the inspiral 
waveform only up to some ^-dependent upper GW fre- 
quency, say oJ^xiy), having the following properties; (ii) 
for v = 0.25, w£tx(0-25) « 0.14 (consistently with [28]) 
which is significantly above the expected range of valid- 
ity of a normal PN approximant, but still significantly 
below the (EOB-estimated 22 ) gravitational wave LSO 



frequency u^so = 0.2114. For v — 0.16, the upper 
bound increases to w^ x (0.16) w 0.17, i.e. around the 
corresponding LSO frequency. For intermediate values 
0.16 < v < 0.25, the situation smoothly interpolates be- 
tween these two cases; (iii) by contrast, as v gets smaller 
than about 0.16, oj^ x (u) will decrease down to values 
of order of 0.05, which are typical of the expected upper 
frequency of validity for a normal PN approximant 23 ; 
(iv) in all cases the range of validity of T4 is limited to 
the inspiral and, contrary to the EOB, does not include 
the plunge; (v) in all cases, T4 exhibits a blow-up of the 
frequency at a finite time. However, this blow up is not 
always the main reason limiting the validity of the ap- 
proximant. For instance, this is the case when v « 0.25, 
but not when v < 0.16. Let us finally emphasize that 
the "enhancement" in the domain of validity of T4 when 
0.16 < v < 0.25 with respect to the normal expected PN 
validity is due to a lucky compensation (which does not 
take place when v < 0.16) between two effects going in 
opposite directions: on the one hand, the bad conver- 
gence of the adiabatic PN expansion; on the other hand, 
the fact that the T4 approximant does not take into ac- 
count nonadiabatic effects (which are quite significant as 
emphasized in [2] and displayed in Fig. 2 of [21]). Our 
present result clarifies the theoretical underpinnings of 
the result found in [65], namely that "deformation" of the 
T4 approximant by spin effects removes the accidental 
nice agreement between T4 and NR. Indeed, one should 
not expect such a chance compensation to be stable un- 
der any deformation of the underlying physics (such as 
additional spins or a varying mass ratio) . 
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For the reasons discussed above, our EOB estimates here use 
as = 25. 

Indeed, Ref. [76] estimated the (3 PN accurate) "PN failure 
point" to be around C 1 / 3 »J 0.3 which corresponds to W22 = 
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